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Since the early days of QCD, it has been argued that inside hadrons quarks organise into sub- 
structures. Among those, popular pictures of various phenomena (including exotica and colour 
superconductivity) give a prominent role to diquarks, i.e. to colour antitriplet combinations of 
two quarks. Using a gauge-invariant setup, which combines the diquark with a static quark, we 
study spatial correlations of the two light quarks inside the diquark. After illustrating the setup, 
we discuss our first results for both the scalar ("good") and the spin one ("bad") diquark channels. 
In a regime in which m K ~ 830 MeV, our data show unambiguously that a scalar diquark forms 
with a size of ~ 0.9 fm. For the vector diquark we are able to put a lower bound of 4. 1(7) fm on 
its size. We also investigate the mass splitting between the good and the bad diquark; our findings 
are compatible with phenomenologically inspired predictions for this quantity. 
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1. Introduction 

The static quark model has provided a rather accurate description of hadronic states over the 
last four decades. This is a rather simple model, since it assumes that quarks inside hadrons are 
uncorrected. Correlations among quarks could however explain remaining features of the spectrum 
that look incompatible with the simple-minded static quark model. In fact a natural possibility is 
that two quarks combine into a colour antitriplet representation. This correlation is energetically 
favourable, since at large distances the resulting chromoelectric field is simply that of an antiquark. 
Colour antitriplet combinations of two quarks are called diquarks. 

Diquarks have been discussed since the early days of QCD (see Ref. [||] for a review). Interest 
in this topic has been recently revived by the increased attention devoted to the high density phase 
of QCD |Q] and to the possible existence of exotic particles [||, |J], since in both cases diquarks 
could play an important role. Diquark-inspired phenomenological models usually start from the 
assumption that those bound states exist. However this assumption needs to be justified from first 
principles. As usual the lattice is an excellent framework for such ab initio calculations. By using 
a gauge-invariant formalism that enables us to study their structure, we show unambiguously that 
diquark states do form and present preliminary results for the mass splittings between the scalar 
and vector diquark channels. Our calculation is done in the quenched theory using Wilson fermions 
and quark masses that correspond to pion mass in the range 830 — 1000 MeV. 

2. Probing the structure of a diquark 

Diquark operators have been discussed in [§]. For convenience, the quantum numbers and the 
colour and flavour representations of the simplest diquark operators are reported in Table. [T| In this 
work we focus on the qcYoYstf (the scalar diquark) and the c/cfq (the vector diquark) combinations. 
Effective colour-spin hamiltonian arguments suggest that the spin zero channel is enhanced with 
respect to the spin one combinations [j|]. For this reason the scalar and the vector diquarks are also 
referred to respectively as the "good" and the "bad" diquarks. 

We now want to investigate whether a di- 
quark bound state forms. The issue is compli- 
cated by the fact that diquarks are coloured ob- 
jects. In order to formulate the problem in a gauge- 
invariant way, we combine the diquark with a 
static quark in a colour singlet representation. Far 
Table 1: The structure of the simplest diquark op- from the source , the influence of the static chro- 
erators (q c is the conjugated spinor of q). m0 electric field becomes less important and the 

physics of the diquark emerges. The static-light 
baryon obtained by combining a static colour source with a light diquark can be studied with the 
density-density correlator method which provides a gauge-invariant technique for studying 
quark distributions inside a hadron. This method has been recently used in Ref. [^, ||] for studying 
the structure of various hadronic states. 

Our setup is sketched in Fig. |l|. At Euclidean time t = we create a static-light baryon in which 
the two light quarks are combined into a candidate diquark channel (either scalar or vector) and we 
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Figure 1: Setup for the computation of the density-density correlators. 



let this state propagate to time t = T /2 (T is the extension of the lattice in the temporal direction), 
where it is annihilated. Far from the source and the sink, at a fixed Euclidean time t (t = T/4 in 
Fig. 1) where the hadron has relaxed to its ground state, we study the density-density correlator of 
the two light quarks 



where the subscripts u and d refer respectively to the up and down quarks, p q = : qyoq '■ and r u 
and f t i are the lattice vector coordinates of the light quarks with respect to the static source. To 
make the analysis of our results clearer, we consider correlators for \r u \ = [r^| = r: the two light 
quarks are on the same spherical shell at distance r from the static source. Since the system is 
spherically symmetric, C only depends on r and the angle 6 = arccos(r„ • fj). If a diquark state is 
favored by QCD dynamics, one would expect that C(6,r) is enhanced at small values of 6, while 
if the quarks are uncorrected the distribution is flat in cos 6. On a lattice the angular distribution is 
distorted by discretisation artifacts, especially for small values of r. To suppress those effects we 
have normalized our distributions to that obtained from uncorrected objects on the same lattice. 

We have measured C(d,r) on 200 quenched configurations at j8 = 6.0 taken from the NERSC 
archive 1 using two degenerate flavours of Wilson fermions at K = 0.153. Those lattice parameters 
correspond to m n ~ 830 MeV and to a lattice spacing a ~ 0.097 fermi. The results are obtained by 
binning in cos 6 within shells of "thickness" 8r = O.lr (so that |r„|/|r^| = 1 ±0.1). The magnitude 
of the correlator decays exponentially with the distance r from the static quark, with a decay length 
of 0.90(2) a (corresponding to a mass of 2.25(5) GeV). Despite this, the signal never dies into 
noise and we are able to study the density-density correlator up to the maximal allowed distance 
r = 8 a on our 16 3 x 32 lattice. In Fig. ^| we report the angular distribution of C(0,r) for a distance 
r = 5. 1 a from the static quark for the scalar diquark and the vector diquark with the z component 
of the spin /- = 1 . An enhancement of the distribution at cos 6 = is evident for the scalar diquark, 

'See http://qcd.nersc.gov. 
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Figure 2: C(0,r) versus cos at distance r = 5.1 a from the static source. 



while the distribution looks flatter for the vector diquark. A study of C(r,d) at larger distances 
r shows the same features. This confirms the expectation that the color field generated by the 
static source does not affect the physics of the diquark once r is large enough, and shows that this 
asymptotic regime can be reached in our calculations. Hence the natural conclusion of this study is 
that a scalar diquark forms, while a vector diquark state, if it exists, is much broader. The density- 
density correlator technique allows also to give a definition of the size of the diquark. At fixed r we 
can look at the correlator as a function of the u — d separation R ui i = 2rsin(0/2). An exponential 
behaviour C(r,R ll( j) oc g-^W^ as seen m ordinary hadrons [|7|], gives a gauge-invariant definition 
of the hadron size Rq. For the scalar diquark, our data are correctly described by an exponential 
ansatz over a wide range of r (see Fig. || for an example). Results for Rq as a function of r are 
reported in Fig. ||. At small separations r the value of Rq depends on the radius r of the shell. The 
extracted values of Rq converge for r > 4 a. Once again, this confirms that the regime in which the 
physics of the diquark system is not disturbed by the field of the static source can be reached on 
our lattice. A safe estimate for the scalar diquark size is Rq = 0.92(9) fm. For the vector channel, 
the fitted size is noisier and seems to increase with r. The value at the maximal possible distance 
r = 8 a can be used to put a lower bound of 4.1(7) fm on the size of this diquark. 

3. Diquark masses 

From colour-spin effective hamiltonian arguments one can estimate the mass of the scalar 
diquark and the mass splitting AM = M vec tor — M sca i ar between the vector and the scalar diquark, 
obtaining respectively M ~ 320 MeV and AM ~ 200 MeV [gp. Since we are using heavy quarks we 
expect these numbers to be different in our case. Again the colour-spin effective hamiltonian gives 
us a handle on what we should expect: if the mass of the degenerate constituent quarks increases, 
the mass M of the scalar diquark increases and the mass splitting AM goes down as 1/M 2 . Within 
our static-light baryon setup the mass of the diquark can not be accessed, since the mass of a 
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Figure 3: The exponential decay of C(r = Figure 4: Fitted diquark size R (in units of a) 
6.15a, R M j) at fixed r for the scalar diquark (the as a function of the distance rja of the diquark 
fit is performed from R ud = 2 a to 8 a). from the static quark. 



j3 


K 


m n (MeV) 


AM (MeV) 


6.0 


0.1530 


830 


115(20) 


5.8 


0.1575 


830 


100(15) 


5.8 


0.1530 


1000 


67(7) 



Table 2: Mass splittings for the scalar and the vec 
tor diquarks at the j3 and K shown. 



hadron with a static quark is ultraviolet divergent. However, the mass difference can be extracted 
by looking at the temporal behaviour of the ratio of propagators of static-light baryons with the 
light quarks combined into the channels we are interested in. 

We study these ratios on a 16 3 x 32 lattice at 
a given pion mass of about 830 MeV at j8 = 6.0 
taking fc = 0.153 and /3 = 5.8 with fc = 0.1575 
to check scaling. At j8 = 5.8 we also calculate 
AM at a larger pion mass of about 1 GeV taking 
K = 0. 153 to check the dependence of AM on the 
quark mass. We use 200 NERSC configurations 
for each value of j3 and K. To reduce the noise 
caused by the static source we used hypercubic smearing [|| [Kj. Our results are reported in 
Table |2] where the errors given are only statistical. Changing the fitting range we obtain values for 
AM that are consistent within these large statistical errors. The following features emerge: (a) at 
fixed m n , the value of the mass splitting is independent of j8, and this suggests that our calculation 
is reasonably free of lattice artifacts; (b) the mass splitting at m n = 830 MeV is smaller than the 
phenomenological estimate; (c) as m n is increased, the mass splitting decreases, (b) and (c) are 
compatible with qualitative expectations from colour-spin effective hamiltonian arguments above. 
Moreover, our results for the mass splittings are also compatible with a similar calculation reported 
at this conference [|TI|], with estimates coming from an instanton liquid model [12] and with an 



earlier lattice calculation in Landau gauge [13] (see also [14]) 



4. Conclusions 



Density-density correlators in the background of a static source provide a gauge-invariant 
method for the study of diquark structure. Within this framework, we study the scalar and the 
vector diquark channels. Our results show that a diquark state forms in the scalar channel, with a 
size of 0.92(9) fm (for a pion mass around 830 MeV). For the vector diquark we can provide a lower 
bound of 4. 1(7) fm on its size. This is the first gauge invariant calculation from first principles that 
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gives evidence for the existence of diquarks as colour antitriplet bound states of two quarks. The 
mass splitting between the scalar and vector diquarks is calculated; our results are compatible with 
extrapolations of phenomenological estimates towards large masses of the constituent quarks. Cur- 
rently we are simulating at a third value of /3 to check scaling and be able to study heavier quarks. 
At the same time we are increasing the statistics at j8 = 6.0 and 5.8 by a factor of 10. By using three 
values of quark masses we should be able to check the 1 /M 2 scaling of the mass splitting between 
scalar and vector diquarks expected at heavy quark masses. Finally we are extending the study 
of the density-density correlation functions and of the mass splitting to the other diquark channels 
listed in Table 1. 
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